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ABSTRACT 

Tuv>dmens>onAl  simulations  of  dun-break  flows 
using  finite  volume  nttthted  and  approximate  Kicmaiui 
solvers  for  computing  the  inter  cell  fluxes  have  drawn 
growing  j  He  nines  because  of  tbnr  robustness  and 
abilities  to  handle  mixed  flows  aiui  discontinuities.  Such 
models  usually  require  complicated  algorithms  for 
treating  source  terms  and  seccrxT  order  w hemes  for 
confuting  the  interccll  flaxes  in  order  to  gun 
numeric  all)  balanced  solutions  jnd  accuracy.  winch 
often  results  in  an  excessively  long  computational  time 
With  a  view  of  developing  an  accurate  and  cfflocut 
model  fee  real-life  applications,  this  paper  proposed  a 
finite  volume  method,  whxrh  uses  the  first-order  KLL 
approximate  Kiemaiui  solver  for  computing  intcrcell 
ll uves  and  adopts  the  coctoteivativc  term  of  the 
momentum  cquMinns  with  one  scarce  term  representing, 
the  driving  forces  in  each  ccputica  Such  treatment  can 
easily  eliminate  numerical  imbalance  between  scarce 
and  lluv  terms  without  introducing  compile  Med 
algorithms.  The  accuracy  and  improvement  in 
cotrpuiMiiinal  efficictK)  of  the  newly  developed  rnxki 
are  deuxostrated  by  iijcorcs  of  a  real-life  tesa  example. 

1.  INTRODUCTION 

The  caustmphic  failure  of  a  dam  often  causes 
ukJcprefei  downstream  Ax»ding.  which  may  directly 
affect  nihility.  deployment,  and  safety  of  the  ami).  A 
2D  dam-hreak  n**A:l  can  capture  both  spatial  and 
temporal  cvoluticci  of  a  potential  dim- break  flood  event 
and  provide  sufficient  details  cs  the  flood,  such  as  flood 
defrhs.  flow  velocities,  aixi  timing  of  the  flood  arnval 
and  rcccssicn  at  specific  local vcis  ard  times.  Svx'h  pieces 
of  information  ar e  crucial  for  military  planning  in  the 
areas  hiving  potential  dmvhreak  risk. 

Numerical  solution  of  two-dinxnsicrul  dam-break 
Dows  has  been  a  great  challenge  to  hydraulic  cngutecrs 
and  researchers,  because  it  often  involves  conplex 
geemetry.  mixed  flow  regimes,  aiui  discontinuities. 
Among  numerous  approaches,  the  use  of  finite  voixme 


methd  w  ith  unstmctured  gyid  and  approximate  Riemann 
solvers  for  computing  th:  intcrcell  fluxes  bus  gained 
ricrcusing  popular  itv.  because  of  its  highly  adaptive 
ability  to  compkex  geometry,  robustness,  abilities  to 
bundle  mixed  (lows  trd  discontinue es  trd  ixAstxvding 
nuss  conservation  property  (c.g..  /ha:  et  al  1996; 
Bmfou  aiui  Ciwcia-Navano  2COU;  Vat  am  ct  al.  2002; 
Yoon  aid  Kang  2004).  However,  when  source  terms  are 
presented  due  to  uneven  buth)n*:iry.  such  methods  may 
create  numerical  imbalance  hccaite  of  the  artificial 
splitting  of  driving  forces  in  the  governing  equations 
between  lluv  ard  source  terms,  which  arc  then  cvaluited 
using  different  methods.  Svx'h  nuiruncal  imhaltttec  can 
bod  an  unph)sical  llow  even  in  a  still  water  test  case,  as 
illustrated  by  Rogers  ct  al  (2Ctt3|  through  a  iwv> 
dimensional  sinailmon  of  a  circular  water  basin  with 
ireven  buthvireuy.  So  far.  many  researchers  have 
attempted  »o  ovcrcctiu  these  probbms  Ter  example, 
Nujk  (1995)  adopted  the  form  of  tbw  governing  equation 
c  which  the  hydrostatic  pressure  force  term  is  cxtrKtcd 
from  the  flux.  Such  treatment  makes  it  possible  to 
discretize  iwo  source  terms  useig  the  san*:  n*Mhcd  and 
drns  to  satisfy  th:  numerical  bulance.  Zhou  ct  al.  (2001 1 
proposed  th:  surface  gradient  method  (SfiM)  for  die 
treatment  of  the  source  terms.  In  tb*  SOM,  water  depth 
at  left  and  right  of  die  interlaces  is  evaluated  bused  un 
die  linear  rcccrstrurtion  of  water  surface  level.  This 
metbted  can  cliiiuxutc  mmcncal  imbalance  problem 
without  introducing  a  coaplcx  algorithm  fer  die  source 
terms.  Rogers  ct  al  (2003 1  used  an  algebraic  approach  in 
which  thw  governing  hyperbolic  system  of  ccettcrvacion 
iiws  is  reformulated  in  terms  of  deviations  away  from  an 
enforced  but  separately  specified  equilibrium  stale  and 
die  numerical  balMVcing  is  achieved  by  uKoeporatng  (he 
resulting  extra  physical  infcnixmcn. 

Th:  use  of  coaplicuted  algorithms  fee  treating  source 
terms  and  the  second-or&r  scheiws  for  computing  th: 
ntcrccil  duxes  often  result  in  an  excessively  long 
compuuncnal  time  (e.g.  Valiam  et  al.  2002;  Yoon  and 
Kang  2004).  In  the  case  of  a  large  scale  problem.  thte  cost 
of  computational  tin*:  can  be  prohibitive.  With  a  view  of 
d:v  eloping  an  Kcuratc  and  cffioerx  model  for  real-life 
applications.  Ying  ct  al  (2006)  presetted  a  fuute  volun*: 
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method.  which  uses  th:  first-order  HLL  uppravimMe 
Kicmann  solver  for  computing  intcnrcll  fluxes  ard  adopts 
the  conservative  form  ot’  the  momentum  equations  with 
one  source  term  representing  th:  driving  forces  in  each 
equiticci  Such  treat  me//.  can  easily  eliminate  iiunuriul 
imbalance  between  source  and  flux  tmrs  without 
introducing  complicated  algorithms.  and  thas  increase 
corrpatMionrd  efficiency.  Nunxncal  tests  have  shmu 
that  tlus  method  is  able  to  sutisfactcnly  predict  oblique 
hvdrautc  jump  r«d  partial  dun- break  flow  <  Ying  et  al. 
2<KK.)  In  this  paper  a  rcul-lifc  test  exanple  with 
coxrpIcA  geoaxtry  (Malpasset  dam-break  case)  is  u*:d  to 
h*ih:r  demonstrate  the  accuracy  and  inprovement  in 
corrpuiMional  efficiency  of  tlus  rew  l>  proposed  method. 


In  the  above  toms  of  dve  Saint  Veiunt  equaiiccK.  th: 
An  Ag  forces  are  represented  by  only  one  term  with  the 
water  surface  gradient  which  makes  it  very  nice  foe 
Beating  the  source  term  because:  1 1  >  the  variation  of 
water  varfrcc  is  generally  much  smoother  than  waer 
d:pth  ard  bottom;  |2)  it  chminrtcs  nimcrcal  mbularcc 
dial  arises  due  to  using  different  methods  to  cvaiuitc 
dnvAg  forces  thit  are  split  between  the  rtuv  and  the 
source  leimt.  As  a  matter  of  fact.  Nujic  <IW5|  proposed 
a  balancing  technique  in  wh>ch  the  pressure  term  dx  to 
water  d:pth  is  extracted  from  the  flux  and  combcrx'd  Ato 
dve  bottom  slcpe  term,  which  actually  resulted  in  th: 
5*1  ux*  form  of  the  governing  cc^utcn  as  Kq.  ( I ). 


2.  COX  KRMNC*  EQUATION’S 

The  twodiirunsiorul  shallow  water  equations  arc 
obtained  by  irxcgrating  dve  Navier-Mokcs  ccpamons  over 
the  flow  &pch  bused  on  several  assumptions  such  as 
hydrc&tMic  pressure  distrihuticA  arcl  sirull  bottom  slope 
The  equaticett  in  corners  aiKin  wd  vector  form  are 
written  as 
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in  which.  U,  MUX  C(U)  and  S(U)  arc  respectively  the 
vectors  of  conserved  variables.  fluxes  in  th:  t  arcl  y 
JirecticA.  and  sources,  defined  as  follows 
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where  A  •  water  depth:  *  and  v  -  velocity  component  in 
the  *  and  >  direction,  respectively;  #  gravitaticeul 
acceleration;  /  “  water  level;  u  -  Manning's  coefficient. 

Fee  convenience.  Fq.  <  1 1  is  often  rewritten  as 
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rt 

where  F-E/  »G  / 


i.  M  MKKK  Al.  ME  I IIOI) 


Fig.  1  Diagrun  fee  control  volum:  definiticn 

Th:  #>vermng  equations  are  doxreti/ed  according  to 
dve  cell  centered  finite  volume  nxthod  on  a  trungulir 
gnd.  as  shown  a  Fig.  I.  Th:  conserved  variables  are 
defined  at  th:  cell  centers  and  represent  the  average 
value  over  each  eelL.  while  the  fluxes  arc  calculated  at 
dve  interfaces  between  cells. 

Iiucgraiing  Eq.  |2|  over  the  /*  cell.  obtains 
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where  O  “  tlu;  area  of  the  f-  cell  Applying  Green’s 
theorem  to  the  second  term  in  Eq.  (5)  yields 
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where  r.  -  boundary  of  the  f  cell  and  ■  -  the  vAit 
coward  vcctcc  rccuul  to  the  boundary.  The  second  term 
in  Eq.  <4 1  can  be  rffreximaicd  as 
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ubeic  the  iul»tnpu  i  and  /  liciklr  lire  f  ceil  x*A  the /* 
edge  ill  the  cell.  rcvpccuvely*.  AT  -  the  length  of  an  edge 
of  a  triangular  cell. 

fherefirc.  lq  <4 » cm  In:  written  as 


Water  surface  gradients  in  the  sourer  teim  an: 
cv  abused  by 
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Aaxtrdmg  to  Godunov  <l95*k  the  varubfc*  oic 
appro  Ainu  ted  as  constant  states  wittun  each  tell  oral  then 
the  llutct  at  innriicet  sec  calculated  by  tolvmg  resultant 
Itiemann  (vublems  that  exte  at  interface*.  In  the  proem 
model.  the  I  ILL  appaiuiuate  Riiminn  solver,  proposed 
by  lianen,  Ln\  and  van  Leer  (llreten  et  al.  I*>343 1.  ts  used 
tu  caLulate  tlte  incrrell  llux.  because  of  its  rubustres* 
and  eo> c  to  urpknctu  According.  to  tin:  I  ILL  Riumuui 
solver.  the  eitavell  tluv  is  given  by 
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where  Ut  and  L*  are  atttrvtd  varuMet  of  the  left  and 
right  sates,  respectively.  aid  .Y»  ted  .V,  an?  left  aid  right 
traveling  wave  speeds,  respectively.  which  are  estimated 
ic  lording  to  the  tollowing  ei|UAtinn*  if-nwcontlkv  and 
Tt.ro  19*51. 


It  is  easy  to  we  Unit  the  resulting  scheme  does  not 
CM&e  nuucricalh  generated  flaw  or  numerical 
iiuhalincc  prcfclan  Lei  us  cousufcr  an  upeii  funnel 
tilled  with  water  ee  rest  Ofevnwisly.  die  water  will  remain 
at  rest  if  nu  dntujharvc  ts  jpptcd  to  die  iLcnun  and 
tc^eniirics.  II  we  apply  above  vclieme  to  smmJjie  this 
erne,  we  can  tee  die  source  term.,  including  water 
surfjee  gradient  term  and  Iriciiuu  tain  and  tlux  lam  in 
lap  <4 1  are  exactly  ccfail  tu  /cm  at  initial  time,  no  nutser 
bow  hid  e levatKin  change*  As  a  result.  die  solutum*  ivf 
/  «d  (J  at  succeeding  turn:  will  remain  die  same  as 
initial  state,  m  other  words,  no  Ikm  is  ewxnericxlly 
generated. 


Fig.  2  I  lipoerapliy  and  locations  uf  measuring  points  fur 
Malpavtci  dam-hreuk  case 


ubeic  l*i  ard  It  are  ruspe.tisvty  the  iducity 
ciuiponenbi  of  the  led  and  the  right  states  in  the  unit 
vector  n  direction:  A,  and  A,  or  c  water  depth  of  the  left 
and  the  ngtu  *utcv  respectively. 
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Stile  cut  for  a  dry  bed  problem  the  wave  spcivh  Si 
and  5k  arc  cstinxiied  iiuvcding  to  the  following 
expression* 
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4.  MODEL  \  AlJDATIONS  USING  AtALFASSE  I 

DAMHKI.XK  CASE 

4.1  Malpassct  Uam- Break  Event 

Itc  Malpasset  dam  was  kxalcd  in  a  current  gorge  ivf 
ftve  Reyran  nver  valley  m  France.  Iz  wjs  a  ut  $  m  high 
arch  dim  with  a  crest  liagth  of  22J  m  and  tic  nuxunum 
reservoir  capicity  of  55*10*  m  In  die  mvnrdute 
ilownurcam  id  tic  dam.  the  Reyran  river  valley  ts  very 
narrow  and  luis  two  consecutive  sharp  tends.  I  lieu  tic 
valley  tVKlem  as  n  goes  downstream  and  eventually 
reaches  the  Hat  plain  <«c  Fig.  2*  IV*  dam  failed  in  I95M 
following  an  excepoiinally  heavy  rain.  After  the  dam 
failure  a  field  survey  was  pcrlomind  to  obtain  tic 
maximum  uuter  levels  along  tfc  Reyran  river  valley  In 
addition,  a  f^v'sical  n»id  wuh  a  scale  of  I  >400  was  iMuh 
to  study  the  dam -br  eak  tluw  in  19<4  IH:  itxixuiaim 
water  level  and  die  timid  wave  armal  tune  at  4  pouits 
alimg  tic  river  valley  were  maisuroi  <sec  lig2| 
Bcvir.ee  of  as  con'pkx  topugi^o  and  availibility  uf 
ircusiircd  data,  the  Xlalfusset  dam-break  case  was 
seteivd  m  a  bcnjluixirk  test  cxairplc  fee  dim- break 
nxdcis  in  the  CADAM  f^cct  (Goutal  IWf.  More 
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detailed  descriptions  about  the  Malpassct  dam-break  test 
cate  cm  be  round  in  literature  (eg.  Goatal  IW, 
Hervouct  2000,  Md  Valiant  cl  al.  2002). 

4.2  Numerical  Simulation* 

In  die  computations^  two  distinct  tnnigular  meshes 
were  respectively  used  in  order  to  investigate  the  impart 
of  mesh  size  co  the  nxdrling  results.  Mesh  A.  as  stown 
in  Fig J.  is  ccetposcd  of  26000  celh  md  highly  refined 
at  the  dim  site,  the  inurcdutcly  dowtwtream  valley,  md 
along  the  downstream  river,  which  allows  the  nxd:l 
accurately  capture  the  details  on  Hooding.  Mesh  B.  os 
shown  in  Fig,4.  is  ccmposcd  of  SK55  cells  aid  not 
excessively  refined  in  ord;r  to  increase  counpuuticcul 
eOlcXicy.  It  should  be  noted  thrt  Mesh  A  was  directly 
odefeed  from  the  CADAM  project  Md  Mesh  B  was 
gewrated  by  tto  commercial  softwMc  SMS.  The  initial 
water  level  in  the  reservoir  was  set  to  he  100  m  above 
sea  level.  The  re*»  of  the  computational  dxnain  was 
considered  as  dry  ted.  The  previous  study  has  stouu 
that  tto  initial  downstream  nver  flow  k  negligible 
becMase  of  its  relatively  auH  flow  rote  comparing  to  the 
flow  earned  hy  the  dim  failure  (Hervouct  2000).  The 
MM»mg‘s  ccctTTcietu  was  set  to  be  0.025  m  ‘  *s  over  the 
entire  con^utational  drenoin  The  time  interval  Jr  - 
(1X125  s  fee  Mesh  Audit  -  0.06  *  for  Mesh  B  The 
calculations  with  a  final  time  r  ”  3000s  require  1  hour  5 
minutes  fee  the  Mesh  A  md  7.6  minutes  foe  the  Me^h  B 
on  a  PC  w  ith  AMD  24  GHz  CPU. 

4J  Result  and  Analysts 

Figs>5  and  6  show  the  fl*xd:d  Mca  and  water  d;pth. 
calculated  using  Mesh  A  and  Mesh  B.  respectively.  It  is 
observed  ttot  the  overall  fleedrd  area  and  water  drpth 
obtained  from  Mesh  A  and  Mesh  B  arc  very  similar. 
deVito  some  tiny  ard  local  differences  in  flooded  area 
md  wrtcr  depth.  It  is  nor  surprising  ttore  Me  such 
difference*.  because  Mesh  B  is  much  coarser  than  Mesh 
A  so  ttot  some  fine  topographic  dMails  such  as 
downstream  river  Me  not  well  represented  in  the 
sinwilotion  using  Mesh  B 

Figs.  7  and  ft  compare  the  ccenputcd  iixixmvum  water 
level  md  wave  front  arrival  tin*:  with  the  measured  dau 
In  these  figures,  good  agreements  between  Mesh  A‘s  ual 
Mesh  B‘s  results  Me  observed  again  Fig.  7  also  stow* 
that  tto  results  of  maximum  water  level  from  the  present 
model  are  very  clove  to  those  frem  tto  model  developed 
hy  Valiant  ct  al.  (2002).  in  which  tto  high  order  schemes 
were  used  fee  bod*  estimating  intcrccll  fluxes  and  time 
advancing.  It  is  import;**  to  note  that  ttoir  colcuLmcoan 
a  quadrangular  mesh  w  ith  10696  cells  reepuret  26  to*** 
on  a  PC  with  Pentium  111  70(1  MHz  CPU.  whereas  the 
present  model  using  Mesh  A  (26000  cells)  c*dy  requires 
4  to***  1 1  nur**es  of  co<i<f**atKinal  time  c*i  a  PC  with 


similar  performance  (Pentium  111  S50  MHz  CPU)  This 
ridicules  tha*.  tto  present  model  cm)  sv4%*ar*ially 
*icrcusc  computaticcul  efficiency,  while  accuracy  is 
preserved. 

5.  CONCLUSIONS 

A  two-dimensional  dom-hreak  model  his  been 
developed  hosed  on  finite  voluix  nxtliod  using 
uiwru:  lured  grids  The  inter  cell  fluxes  are  cvaluited 
based  on  the  HU.  appeaxinxite  Rienxinn  solver.  Tto 
nomcr/.um  equation  used  in  the  model  has  only  one 
source  term  representing  the  driving  forces.  This 
approach  can  successfully  eliminate  tto  numerical 
imbalance  between  source  and  tlux  tarns.  Tto  model  is 
validated  against  lahoratccy  data  from  a  real-life  dim- 
break  case.  The  results  show  thit  the  simulations  using 
cittor  MeA  A  or  Mdt  B  arc  c rouble  of  correctly 
capturing  major  hydrodynamic  behaviors  of  the  flood 
event,  such  as  nxixinwsn  wr*er  level  and  wave  front 
arnvol  time.  The  comparison  of  ccenputMional  times 
with  other  models  tv*ttor  demonstrate*  that  the  present 
nttdei  is  more  cfficia*  thin  conventiccul  unstructured 
nxdebs  using  high-order  schemes  for  cstimaiing  fluxes 
and  ccmplicated  upwind  algorithms  for  treainx*nt  of 
sc**ve  term.  Such  an  unstru:  lured  nixkT  allows  using 
arhitrMy  cecnpwtation.il  domain  and  meshes  with  local 
rctineiixnt.  This  provides  the  more  efficient  way  to 
otoain  local  fine  details  on  flooding  than  stnicturcd 
iwdek 

To  sum  up.  tto  proposed  model  is  able  to  correctly 
ensure  both  spatial  m* 3  temporal  evolution  of  a  potential 
dam-break  flood  event  and  provide  crucial  information 
for  military  planning  in  tto  areas  having  pkxeufial  dim- 
bccok  risk. 
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